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The strong isospin-breaking correction, fi st , which appears in estimates of the Standard Model 
f*^ ■ value for the direct CP- violating ratio e'/e, is evaluated to next-to-leading order (NLO) in the chiral 

expansion using Chiral Perturbation Theory. The relevant linear combinations of the unknown NLO 
CP-odd weak low-energy constants (LEC's) which, in combination with 1-loop and strong LEC con- 
tributions, are required for a complete determination at this order, are estimated using two different 
models. It is found that, to NLO, fl 3 t = 0.08 ± 0.05, significantly reduced from the "standard" 
value, 0.25 ± 0.08, employed in recent analyses. The potentially significant numerical impact of this 
decrease on Standard Model predictions for e'/e, associated with the decreased cancellation between 
C****" , gluonic penguin and electroweak penguin contributions, is also discussed. 
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^ ; I. INTRODUCTION 

The recent improved experimental results for the ratio of direct to indirect CP-violation parameters, e'/e, obtained 
by both the KTEV and NA48 collaborations [1,2], have spurred on continuing efforts to reduce the theoretical 
uncertainties in our expectations for the value of e'/e in the Standard Model. While short-distance effects are under 
r"| • control (the Wilson coefficients of the effective weak Hamiltonian being known to two-loop order) [3] , there remains 
ij significant uncertainty in the theoretical calculation of the long-distance K — > 7T7t hadronic matrix elements. These 

have been estimated using a number of techniques and models in Refs. [4-10] (for recent reviews, see Refs. [11,12]). 
" An important ingredient in the calculation of the hadronic matrix elements is the inclusion of isospin-breaking 
(IB) effects. Strong isospin breaking arising from the up/down quark mass difference, 5m = nid — m u , induces a 
AI = 3/2 contribution to the K — » irir matrix elements of the gluonic penguin operator, Qq 1 , which, in the isospin 
limit, is pure AI = 1/2. This 'leakage' of octet (AI = 1/2) strength into the AI = 3/2 component of the K° — ► tttt 
decay amplitudes has the effect of reducing the magnitude of the Qq contribution to e'/e which one obtains in the 
isospin-conserving (IC) limit. This is conventionally represented by a multiplicative factor 1 — Q st applied to the 
AI —1/2 matrix element. Explicitly, one writes [3] 
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with 



p(l/2) _ p(3/2) 



(1) 



p(i/2) = ^ yi (Qi) (i - n st ), 

p(3/ 2 ) = I^ yi(Qi ) 2; (2) 

UJ — 

where lo ~ 1/22, the yi are the parts of the Wilson coefficients associated with the top quark (and hence, through the 
corresponding CKM matrix elements, with direct CP-violation), and the subscripts and 2 denote the isospin of the 
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7T7T final state. Note that we have dropped overall factors in Eq. 1 in order to highlight the dependence on Q st . In 
typical analyses one finds a significant cancellation between the AI =1/2 and AI = 3/2 contributions, and therefore 
a non-trivial sensitivity to Q s t ■ 

The IB correction, Q st , is obtained as follows. Writing the isospin decomposition of the K° — ► tttt decay amplitudes 

as 



A + . = ^ l -A e i *° + ±=A 2 e i *\ (3) 

where A$ and A 2 are, respectively, the (in general complex- valued) AI = 1/2 and AI = 3/2 amplitudes, and fa = Sf* 
are the usual ir — it scattering phases, fl s t is given by the ratio 

^ = u^A Q (4) 

where SA 2 is the octet leakage contribution to the AI = 3/2 amplitude, and u> = (ReA 2 /~ReAo) — 1/22.2 (reflecting 
the AI =1/2 rule enhancement). 2 

The possibility that the octet leakage contribution to the matrix element of Qq induced by 8m ^ and electric 
charge differences could have a significant impact on estimates of e'/e was first discussed in Refs. [13,14]. At leading 
order (LO) in the chiral expansion, the leakage contribution to A 2 is saturated by ir°-r) mixing and the kinematic effect 
produced by the K°-K ± mass splitting and the momentum-dependence of the LO weak vertices. This leads to the 
well-known LO value Q st = 0.13. The difference between this LO value and the conventional value employed in recent 
analyses of e'/e, fi st = 0.25 ± 0.08 [14], results from an estimate of those NLO effects mediated by the rf through the 
K — ► ttt]' transition and 77-77' mixing, which effects would be expected to be dominant in the large-7V c limit. In the 
framework of the conventional low-energy effective theory employed in this paper (that involving only the tt, K and r\ 
degrees of freedom, in which the rf and other higher resonances have been integrated out), such effects correspond to 
contributions to the NLO weak LEC's (to be discussed below), although use of phenomenological values for the octet- 
singlet mixing angle also effectively incorporates NLO //-mediated contributions proportional to the renormaliscd 
strong LEC, Ly, of Gasser and Leutwyler [15]. As has been recently pointed out, however, other NLO contributions 
might also be important. An example is that discussed in Ref. [16]. If one considers the effect of NLO strong dressing 
on the external legs in K — > tttt, there is a large NLO ^'-induced contribution (proportional to L T 7 ) associated with 
treating 7r°-?7 mixing at NLO. This mixing contribution, however, always occurs in the fixed combination 3L T 7 + L'^ 
(see, for example, the expressions for the angles, Q\ and 6 2 , describing NLO mixing, given in Ref. [17]). As pointed 
out in Ref. [16], there is a strong numerical cancellation between the ^'-induced Uj and the scalar-resonance-induced 
ig contributions, clearly demonstrating the importance of including NLO contributions other than those induced by 
the rj . The effect of this cancellation was found, in Ref. [16], to lower Q st to 0.16 ± 0.03 (thus increasing the Standard 
Model prediction of e'/e by about 21% [16]). The possibility of additional non-r/-induced NLO contributions to 
£l st was also recently explored in Ref. [18]. The focus of this work was on IB contributions associated with NLO 
weak LEC's, and the numerical results suggest the possibility of very large corrections to Q st associated with scalar 
meson exchange. However, as we discuss below, the results of Ref. [18] suffer from important technical shortcomings 
which make them numerically unreliable. The significance of the possible reduction in the value of f2 s t suggested by 
Rcfs. [16,18] is obvious from Eqs. 1 and 2. 

Chiral perturbation theory (ChPT) provides a natural framework for the calculation of Q st since it ensures that all 
contributions of a given chiral order may be obtained in a computationally straightforward manner. The complete 
set of NLO contributions is a sum of NLO strong LEC, one-loop, and NLO weak LEC contributions, each of which 
is separately renormalisation-scale-dependent, divergent, and therefore unphysical. The sum of these contributions 
is, however, necessarily finite and scale- independent. The cancellation of both divergences and scale-dependence in 
the final result provides a highly non-trivial check of the explicit calculations. The set of Feynmann graphs to be 
evaluated is shown in Fig. 1, where Fig. 1(a) represents the LO, Fig. 1(c) the NLO strong LEC, Figs. 1(b) and 1 (d)-(g) 
the one-loop, and Fig. 1(h) the NLO weak LEC contributions, respectively. In Ref. [16] only the NLO contributions 
associated with external line dressing (Figs. l(b,c)) were considered, while Ref. [18] examined only the contribution 
of Fig. 1(h). 



2 We follow standard phase conventions in which Aq and Ai are real in the absence of CP-violation. 
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In Refs. [19,20] it was shown that, in the CP-even sector, the inclusion of one-loop isospin-breaking corrections in 
K — ► 7T7T decay amplitudes acts to decrease the magnitude of the octet leakage contribution to the AI = 3/2 amplitude. 
This decrease ensures that the AI = 1/2 rule enhancement obtained in an IC analysis is accurate to better than 10%. 
In the present paper we obtain a complete NLO (0(p 2 Sm)) determination of Q st in ChPT, and show that the same 
qualitative situation holds in the CP-odd sector. This means, according to Eq. 4, that one should expect a further 
reduction of fl s t beyond that associated with the NLO strong mixing effects studied in Rcf. [16]. The absence of 
reliable values for the CP-odd weak LEC's is the crucial stumbling block on the way to an accurate numerical result 
for flgt at NLO, the NLO strong LEC and loop corrections being, as we will see below, well-determined. In the present 
work we appeal to two models (the weak deformation model, and the chiral quark model) of the contribution of the 
gluonic penguin operator, Q 6 , to the CP-odd weak LEC's in order to probe the probable scale of the model-dependence 
in our estimates of Q st . It is important to emphasize that, while the estimates of the weak LEC contributions to 
il s t are model-dependent, the one-loop and strong LEC contributions (discussed below) are, though scale-dependent, 
model- independent. As we will see, the combination of these model-independent NLO contributions is, at typical 
hadronic scales, rather large, and negative, suggesting a significant reduction of Q st as compared to the conventional 
value. 

The rest of the paper is organised as follows: In the next section we briefly review the chiral Lagrangian approach 
to the calculation of non-lcptonic kaon decay amplitudes and discuss the models employed for the relevant NLO weak 
LEC combinations. In Section III we present our numerical results for fl st . The impact of our findings on theoretical 
estimates for the value of e'/e in the Standard Model are discussed and conclusions are presented in Section IV. 



II. CP-ODD K° -> 7T7T DECAY AMPLITUDES IN CHIRAL PERTURBATION THEORY 

The diagrams which have to be calculated to obtain the K — > tttt decay amplitudes to NLO are, as noted above, 
those given in Fig. 1. We now briefly review the ingredients needed for these calculations, referring the reader to 
Ref. [19-21] for the technical details. The low-energy representation of the non-leptonic weak interactions is obtained 
from the effective chiral Lagrangian, £\y, which was written to LO in Ref. [22], and up to NLO in Ref. [23] (or, in 
equivalent reduced forms, in Refs. [24,25]). (In the CP-odd case, the gluonic penguin operator, which is the focus of 
the present work, and which, together with the electroweak penguin operator, dominates e'/e in the Standard Model, 
is pure octet; we, therefore, need only the octet components of Cw) We work with a form of the effective weak chiral 
Lagrangian in which the weak mass term appearing at LO [22] has been rotated away [23] . The LO (second order in 
the chiral counting) part of the octet Lagrangian, C W ( 8 y is thus given, in the absence of external fields, by [22] 

£$8) = 4T± [^d^d»U] (5) 

where U = cxp(iA • n/F), with Tr a the usual octet of pseudoscalar fields, and F is the pion decay constant in the 
chiral limit. The superscripts ± label the CP-even and odd cases respectively, with A+ = A 6 and A~ = A7. thus 
represents the LO CP-odd octet weak coupling strength. 

The NLO octet weak effective chiral Lagrangian is similarly given by either [24] 

± 37 

or [23,25] 

48 

where the operators which contribute to the K — ► itit amplitudes correspond to i = {5, 6, ...13} and i = 
{1, ...5, 10, ...15, 32, ...40} in Eqs. 6 and 7 respectively. In quoting our results below, we will employ the notation 
of Eq. 6, and hence work with the operator basis given by 

5 (X ± ) = Tr[\ ± {S,L^}] 
6 (A ± ) =Tir[A ± L M ] Tr[SL»] 
7 (A ± ) = Tr[A ± 5]Tr[i^] 
8 (A ± ) = Tr[A ± J L^]Tr[5] 
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9 (X ± ) = Tr[\ ± {P,L^}] 
O la (\ ± ) = Tr[\ ± S 2 ] 
11 (\ ± ) = Tr[X ± S]Tr[S] 
12 (X ± ) = Tr[\ ± P 2 ] 

13 (A ± ) = Tr[A ± P]Tr[P] (8) 

where = iWd^U, S = x f t/ + W X and P = i( X ^U - W X ), with X = 2B M q (where M q is the quark mass matrix). 
Note that in Eq. 6 the weak LEC's are expressed as products of factors c 2 /F 2 and N~ . The latter will, henceforth, 
be referred to as reduced CP-odd NLO weak LEC's. (As we will see below, this reduced form has certain advantages 
for estimates of fl s t-) 

The remaining ingredient needed in order to calculate the diagrams of Fig. 1 is the strong chiral Lagrangian. We 
use the standard form of Gasscr and Leutwyler given by Cs — + + • • • , where the superscripts indicate the 
chiral order and, in the absence of external fields, one has [15] 

4 2) = X Ir^UPtf] + ^Trlxf/t + U X \ (9) 

4 4) = ^(Tr^t/^C/t]) 2 + L 2 Tr[d t ,Ud u U^} Tr[d" Ud v U^] + L^d^Ud^ d v Ud v U^] 

+L 4 Tr[d M £/^[/ t ] Tr[ X Ui + U X ^\ + L^d^Ud^ { X U^ + U X ^)] + i 6 (Tr[ X C/ t + U X ^]f 
+L 7 (Tr[ X tf _ f/x t ]) 2 + L8 Tr[ x [/t xC/ t + U} A U} A] + H 2 Tr{ XX ^} , (10) 

where {Li}, F and Bq are the usual strong LEC's, in the notation of Ref. [15]. Recall that, when using dimensional 
regularisation, the NLO LEC's, {Li}, are formally divergent and have a Laurent expansion in d — 4 (where d is the 
spacetime dimension) of the form 

L t = r t \ + U i +L ( i - 1) (d-4) + ... (11) 

where 

2 



A 



32tt 2 



d-A 



+ j E - 1 - hi(47r) 



(12) 



The {L\} are the usual scale-dependent renormaliscd versions of the Li [15], for which we employ the values found 
in Ref. [26], while the {F^} are constant coefficients (frequently called scaling coefficients) which govern the scale 
dependence of the L\. The contribute first to physical observables at ncxt-to-next-to-leading (sixth) order 

(NNLO), through one- loop graphs involving a single NLO vertex proportional to Li and, as such, are on a similar 
footing as the LEC's present in . (The NLO weak LEC's, {Nf}, in Eq. 6, of course, have a similar expansion.) 

The formal difference between the K° — > tttt vertices extracted from Eq. 5 using A + and A~ is the switch — > ic^. 
This is also, therefore, the only difference between the corresponding LO CP-even and CP-odd decay amplitudes. 
Since the NLO strong LEC and one-loop contributions to these amplitudes involve a single LO weak vertex, one 
readily sees, from Fig. 1, that the substitution — > ic^ also converts the CP-even version of these contributions into 
the corresponding CP-odd version. The substitution c^N^ — > ic^ N~ , similarly accomplishes the CP-even— > CP-odd 
conversion of the weak LEC contributions, Fig. 1(h). If one considers the ratio of the NLO contributions to the LO 
contributions, therefore, the only difference between CP-even and CP-odd cases for the K° — > tttt amplitudes is the 
difference in the numerical values and physical interpretation of the renormalised reduced weak LEC's, N[. Thus, 
for example, the contributions to ImSA2 arising from the diagrams of Fig. l(b)-(g) are immediately obtained from 
the first row of Table II of Ref. [19] by multiplying the entries by ic^. The formal contributions from the CP-odd 
NLO weak LEC's, [AjJ^lec an d [^^Wlbc, arising from Fig. 1(h) are given in the Appendix in terms of the of 
Eq. 6. Unlike the case of the CP-even sector, where linear combinations of octet NLO weak LEC's corresponding to 
isospin-conserving contributions were fit to the available K — ► tttt and K — ► 7T7r7r data in Ref. [27], both the CP-odd 
isospin-breaking and isospin-conserving NLO weak LEC combinations arc unknown. As there is not sufficient data 
available to perform fits analogous to those in the CP-even sector, it is necessary to resort to models to estimate the 
weak LEC contributions. 

In general the numerical value of the isospin-conserving combination of NLO weak LEC's, which appears in 
[^o]wlbc (Eq. Al of the appendix), can be determined from the expressions for the hadronic matrix elements 
at NLO (for which many calculations exist). However the IB combination which appears in [iJAJvyLBC cannot be 
so determined. In what follows we consider three models which could be used to estimate the NLO weak LEC con- 
tribution to fl s t- These are the weak deformation model (WDM) of Ref. [24] from which direct NLO weak LEC 
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estimates are available, the chiral quark model (%QM) as implemented by the Trieste group in Refs. [8-10] for which 
all the necessary ingredients required to calculate the NLO weak LEC contributions to Im^4o an d hnSA2 are readily 
available, and the "scalar saturation model" of Ref. [18], which combines the factorisation approximation with the 
assumption of scalar meson exchange saturation of the relevant strong LEC's. 

The weak deformation model of Ref. [24] proceeds from the observation that the LO weak chiral Lagrangian of 
Eq. 5 can be generated from the LO strong chiral Lagrangian of Eq. 9 by a simple "topological deformation" . The 
model hypothesis is that this same deformation can be used to generate the entire AS = 1 chiral Lagrangian. The 
WDM thus provides no information about the LO weak LEC values, c ± , but gives explicit expressions for the reduced 
NLO weak LEC's, Nf, in terms of the NLO strong LEC's, L t . For the NLO weak LEC's relevant to K —> rnr we 
have [24] 

3 

[N 5 ]wDM = — 2 ["^6 }WDM = ~[N 7 ]\YDM = ~ L5 

[Nz]wdm = 4L 4 + 2i 5 . (13) 

All other NLO weak LEC's appearing in Eq. 6 vanish in the WDM. Since, in these relations, the divergent parts 
of the {Li} do not generate the correct divergent parts of the {N~}, one must interpret Eqs. 13 as applying to the 
renormalised versions of the LEC's. Moreover, since the scaling of the weak LEC's is not correctly given by that of 
the strong LEC's, Eqs. 13 can be taken to hold only at a single scale, which is assumed to be a typical hadronic scale, 
Hh- Assuming resonance dominance suggests fih ~ m p . Note also that, because of the scale dependence of the weak 
LEC's, even though, in the WDM, the remaining N[ vanish at the assumed matching scale /j,h, they are non-zero 
at other scales. Numerical values for the IC and IB combinations of NLO weak LEC's relevant to K° — > irir in the 
WDM can be found in the last column of Table I, where we have used conventional values for the renormalised strong 
LEC's, L\ and L\ at [i = m p taken from Ref. [26]. 

In another approach, the hadronic matrix elements of the four-quark operators relevant for non-leptonic kaon decay 
have been estimated in the xQM- by the Trieste group [8-10]. In particular, expressions for the gluonic penguin 
contribution to the LO weak LEC's were obtained in the second of Refs. [10], and preliminary estimates of the 
contribution to NLO weak LEC's were given in Ref. [8]. As pointed out by the authors of Ref. [8], however, the latter 
expressions contain errors and arc not to be used [28] . We have used the updated results of Ref. [9] , which give the 
contributions of Q§ to Aq to NLO, to extract the CP-odd weak LEC combination which enters Irm4o- To obtain the IB 
combination of NLO weak LEC's entering SA 2 , we have calculated the matrix element (7r + 7r° | Q 6 1 = (^3/2)5^2 
in the xQM using the formalism described in Ref. [9], together with the corrected results for the basic ingredients 
required to compute such matrix elements in the model given there. The model parameters which enter these estimates 
are the constituent quark mass, M, the vacuum quark condensate, and the gluon condensate. These parameters were 
constrained to reproduce the AI =1/2 Rule in the CP-even sector in Ref. [9] using a matching scale A x = 0.8 GeV 
(the scale at which the scale-dependence of the chiral loops and that of the short-distance expressions is roughly 
matched in the model). The fitted parameters have the values M = 0.2 ± 0.02 GeV, (qq) = (-0.240t?o GeV) 3 , and 
(a s GG/-K) — (0.334 ± 0.004 GeV) 4 . The resulting values for the NLO weak LEC combinations are presented in the 
first column of Table I. Numerical estimates for the NLO weak LEC contributions to Im^o and lmSA2 are presented 
in the next section. 

The third approach to estimating the CP-odd weak LEC's is that of Ref. [18]. In this reference, one begins with 
the "factorisation approximation" for Qq 3 , in which the low-energy representation of Qq is assumed to be given by 
the product of the low-energy representations of the scalar densities of which the unrenormalised operator is formally 
a product. Since the low-energy representation of each such density is obtained by taking the derivative of £s[x, x^], 
with respect to the appropriate component of \ or X 1 * (treated here as external sources) [15], the model version of 



3 The factorisation approximation becomes exact in the limit of large N c . In this limit, taking Qq for example, if one 
renormalises the two densities of which Qq is a formal product, then one will also have renormalised the four-quark operator 
Qq. For two such renormalised densities, J(x) and J'{y), at different points x, y, one can straightforwardly construct the low- 
energy representation of the product J(x)J'(y) using standard methods. In general, of course, the low-energy representation of 
such a product is not simply the product of the low-energy representations of the individual densities, but also contains seagull 
terms. It turns out that, for Qq, the NLO part of this representation, which is the part investigated by the authors of Ref. [18], 
does indeed contains seagulls. Since these diverge as x — * y, it is necessary to interpret the factorisation approximation as 
corresponding to an approximate low-energy representation of Qq obtained by dropping the seagull terms in the low-energy 
representation of lim x -> y J(x)J'(y), i.e., to one obtained by taking simply the product of the low-energy representations of the 
two densities. 
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the low-energy representation of Q§ becomes the product of two such derivatives. This product can be organized by 

(2) 

chiral order. The LO term (second order in chiral counting), arises from the product of the derivatives of £ s and 
£g\ and leads to the conventional factorisation approximation for c% , 

c 2 - = ^V us V ud lmC 6 (16B 2 F 2 L 5 ) (14) 

where Cq is the Wilson coefficient accompanying Q e in the expressions for the effective weak Hamiltonian. 4 At NLO, 
having dropped the seagull contributions associated with the second derivative of with respect to the sources, one 
is left with two types of terms, those involving a product of two derivatives of , and those involving a product of 
one derivative of and one of . Let us denote the resulting contributions to the factorisation approximation for 
the (non-reduced) NLO weak LEC's, E~ , (as employed in Ref. [18]) by and [-Ej~] 2 , respectively. The authors 

of Ref. [18] then employ a model in which the relevant renormalised fourth order and sixth order strong LEC's are 
assumed to be saturated by scalar resonance exchange. There is, however, a problem with the approach of Ref. [18], 
associated with the [E~] contributions. To understand the origin of this problem, consider, for example, the results 
of Eq. 17 of Ref. [18], translated into our notation: 

[ E i] i = [ E s] x = [-£ 5 _ ] x = ^V us V ud ImC 6 (32B 2 )L 2 . (15) 

As noted above, L% has a Laurent expansion of the form given by Eq. 11. Note, first, that this means that the 
contributions [E~] of Eq. 15 begin at 0[l/(d — 4) 2 ], in contrast to the actual , whose Laurent expansions begin 
at 0[l/(d — 4)]. While one might plausibly ignore this discrepancy, arguing that only the finite parts of the expressions 
at some hadronic scale are to be used in any case, a related problem remains, even for the finite parts. Explicitly, the 
fact that both factors of L$ in Eq. 15 contain a l/(d— 4) term means that the finite part of L\ is not [Lg] 2 , as assumed 
in Ref. [18], but rather [Lg(/z))] 2 + 5/(3847r 2 )L 8 ~ 1 ^(^i), where the explicit value ofTg, given in Ref. [15] has been used. 
Since, as explained above, the h\ ^ arc on the same footing as the sixth order strong LEC's which enter [-Ej 7 ^' the 
model expressions of Ref. [18] for the sum [£""] x + 2 are numerically incomplete. The model, moreover, provides 

no means of estimating the L\~ r \ making it impossible to correct this defect. In view of this problem, we conclude 
that, at present, it is not possible to estimate the NLO weak CP-odd LEC's using the factorisation approximation. 

Given the problem just discussed with the numerical estimates of Ref. [18], we restrict our attention to the WDM 
and xQM in obtaining estimates for the NLO weak LEC contributions to ImA and Im<L4 2 . The resulting values 
for Q st will be given in the next section. Since the xQM is a microscopic model, and the WDM is not, we will take 
the value of Q s t obtained using the former model as our central value, and use the deviation from this central value 
of the result obtained using the WDM as a minimal measure of the theoretical uncertainty in our prediction for fl st 
associated with the model dependence of the weak NLO LEC's. 

III. NUMERICAL RESULTS 

The isospin-breaking correction to the gluonic penguin operator, Q 6 , evaluated to 0(p 2 Sm), can be written in terms 
of its LO (O(Sm j) value and NLO corrections as 

n[f [1 + R 2 - Ro] (16) 

with 



)(2 ) _ \/2 B Q 8m 

6w (m 2 K - ml) 



where vtik — {m K o + m K +)/2. (Note that the result of Eq. 17 is unambiguous and independent of the LO weak 
coupling Cj .) The NLO corrections, are given by 



4 In ChPT, c 2 is finite, and scale-independent, whereas L5 is divergent. To make sense of this relation one, therefore, usually 
assumes that L5 is to be replaced by its renormalised value, L^(fj,), evaluated at some typical hadronic scale fj, = Hh- 
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ff _ Im 4 NL °^ Im S4 NL °' ND) riR , 

ImA^ ' 2 " In,<» (8) 

where the superscript (NLO, ND) indicates the sum of non-dispersion NLO contributions (involving NLO weak and 
strong LEC's and the non-dispersive parts of loop graphs). Note that, in Eq. 18, R 2 arises from IB effects, whereas R 
is purely isospin-conserving. (The IB correction to Rq would generate a contribution to Q st of O ([<5t7j] 2 ), and thus is 
beyond the scope of the present work.) To separate the model- independent contributions associated with strong LEC 
and loop effects (Figs. 1(b)- (g)) from those of the model-dependent NLO weak LEC terms (Fig. 1(h)), it is convenient 
to further expand Ri as 

Ri = R {non-WLEC) +R {WLEC) (1Q) 

where the superscripts indicate NLO weak LEC (WLEC) and one-loop-plus-strong-LEC (non-WLEC) contributions 
respectively. We begin our discussion with the non-WLEC contributions. These arc model-independent and unam- 
biguous, albeit rcnormalisation-scale-dcpcndent, since the NLO contributions concerned all involve exactly one weak 
0(p 2 ) vertex. The resulting overall factor of in the non-WLEC part of the numerator of Eq. 18 therefore cancels 
with the corresponding factor in the denominator. This cancellation removes all of the short-distance uncertainties 
(Wilson coefficients, CKM matrix elements, etc.) contained in . In addition, because the LO coupling strength 
cancels and the strong vertices (if any) are identical for the CP-even and CP-odd cases, diagram-by-diagram, the 
non-WLEC contributions to Vt st are 'universal', that is, they are the same for CP-even and CP-odd cases. To evaluate 
these contributions we use as numerical input the values = 135 MeV, % = 495 McV, = 549 MeV, and 

Bo Sm = ( md - m A m l = 5552 ± 674 MeV 2 (20) 



(2) 

(as determined by Leutwyler in Rcf. [29]). With these values, we have the usual result fl^.' = 0.128 w 0.13. The NLO 
non-WLEC contributions to Rq an d R2 are given in Table II. The results are presented at two different renormalisation 
scales, fi — m r] and fi = m p , in order to display explicitly the scale dependence of the non-WLEC contributions. When 
using the WDM estimate for the NLO weak LEC's, we employ [i = m p (consistent with the expectations of resonance 
saturation), and when using the xQM estimates, /x = 0.8 GeV (the matching scale employed in Ref. [9] in obtaining 
fits for the xQM parameters). 

It is immediately apparent that the non-WLEC contributions to Rq and R2 both act to reduce Q s t as compared to 
its LO value. The Ri non-WLEC contribution is weakly scale dependent and, being "universal", follows immediately 
from the corresponding CP-even results of Ref. [19]. Although for scales \i ~ m p R ^ lon - WLEC ) - ls positive, and 
hence acts to lower ft st , the scale dependence, in this case, is significantly stronger. The increase in the magnitude 
of ImAo associated with the loop contributions is what one would expect given the attractive final state interactions 
(FSI) in the I — channel, and is analogous to the Aq FSI enhancement discussed previously for the CP-even 
case [30,31,27]. The effect of FSI on the CP-odd amplitudes has also been recently discussed in Refs. [32-34]. In 
Refs. [32,34] it is argued that the value of lmA obtained from approaches which do not generate the final state irn 
I = 0, 2 phases for the / — 0, 2 K — > tttt amplitudes should be enhanced by FSI by a factor of ~ 1.55, while the 
value of ImA 2 should be suppressed by the weakly repulsive 1 = 2 FSI by a factor of ~ 0.92. The numerical values of 
the enhancement /suppression are obtained using the Omnes representation for the amplitude, and correspond to the 
subtraction point s = 0, for which the ChPT representation of the amplitude is presumed to be accurate. It should be 
borne in mind that the I = 0, 2 FSI, corresponding to Fig. 1(f), are already correctly included in our calculations, up 
to NLO in the chiral expansion. That the L = FSI, for example, produce a significant enhancement of ImAo, follows 
from the known FSI enhancement of Ao in the CP-even case [27] and the "universality" of the one-loop contributions. 
One should also note that the Omnes function part of the representation of the amplitudes does not incorporate 
all of the NLO effects; some NLO effects remain in the polynomial prefactor. Thus the question of interest to us, 
namely whether the complete set of NLO contributions raises or lowers ImA relative to its LO value (i.e., whether 
Rq is positive or negative) is not determined solely by the character of the loop contributions; it is perfectly possible, 
in principle, for the NLO LEC contributions to be sufficiently negative that the full NLO determination of Rq is 
negative, even in the presence of the attractive FSI phases. For the models we have considered, this is not the case, 
and the combination of non-WLEC and WLEC contributions to Ro is positive, leading to a suppression of fi st below 
its LO value. It is important to note that, so long as one adheres to the convention of incorporating the effect of the 
1 = 2 leakage contribution by means of a multiplicative correction factor applied to the contribution to the 1 = 
amplitude, there is an amplification effect at work in the gluonic penguin contribution to e'/e associated with the 
NLO contributions to ImAo: the more NLO effects increase the isospin-conserving contribution to ImAo, the more 
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they simultaneously decrease Q s t- Since the contribution to e'/e, including IB, is proportional to the product of the 
isospin-conserving contribution and the factor 1 — Q st , both effects serve to enhance the Qe contribution to e'/e. 

The NLO weak LEC contributions to Rq and i?2 are estimated using the models described in the previous section. 
The numerical results are displayed in Table III. We note first that the contributions to R I ^ VLEC - 1 and R 2 WLEC - 1 
in the WDM are identical, and hence cancel in the difference, R2 — Ro, entering Eq. 16. In the xQM, the WLEC 
contributions to fl st are positive. Indeed the results of Tables II and III show significant cancellation between the 
WLEC and non-WLEC contributions in the xQM. 

(2) 

The total NLO correction factor, 1 + R2 — Rq which multiplies fl st , resulting from the combination of WLEC and 
non-WLEC contributions, is 

1 + R2 - Ra -\0.27 (WDM). (21) 

Taking the xQM result as a central value, and the deviation of the WDM result from this central value as a minimal 
measure of the mo del- dependence of our result, we find the IB correction to the gluonic penguin contribution to e'/e 
to be 



fl st = 0.08 ±0.05 ±0.01 (22) 

where the first error represents the uncertainty associated with the model dependence of the NLO weak LEC's, and 
the second the uncertainty in the input value of B Sm. The central value in Eq. 22 is significantly lower than both 
the conventionally-employed value, 0.25 ± 0.08, and the result of Ref. [16], 0.16 ± 0.03. That the relative uncertainty 
increases from about 30% to 62% is a reflection of the uncertain state of our knowledge of the NLO weak LEC's. We 
emphasize, however, that, regardless of the actual value of the weak NLO LEC contributions, the model-independent 
one-loop and strong LEC contributions, which are unambiguous, produce a significant reduction of U st for any plausible 
choice of hadronic scale. As such, the inclusion of the loop contributions is crucial to any attempt to evaluate f2 st 
beyond LO. For the models considered for the weak NLO LEC's, the net effect is to drive fl s t significantly below its 
LO value. 



IV. DISCUSSION AND CONCLUSIONS 



As noted above, the significant cancellation between gluonic penguin and electroweak penguin contributions means 
that predictions for the value of e'/e in the Standard Model can depend rather sensitively on fi s t. The exact degree 
of sensitivity, of course, depends on the relative size of these two dominant contributions, on which there is, as of yet, 
no clear theoretical concensus. In order to illustrate the impact of the decrease of fl s t from the conventional central 
value, 0.25, to 0.08 ± 0.05, let us use the rough approximation to Eq. 1 discussed in Ref. [5] 

1 a [B a (l - fi rt ) - 0.4B 8 ] (23) 

(where we have dropped an overall constant multiplicative factor irrelevant to the present discussion). Maintaining 
the constraint, B & > Bs, imposed by the Munich group [5,12], and using the values B 6 = 1.0 ± 0.3 and B s = 0.8 ± 0.2 
employed by them, we find the results shown in Table IV. The range of values for Bq, Bg covered in the Table is the 
same as that in Table 3 of the second of Refs. [12], from which the values for e'/e corresponding to Q st = 0.25 have 
also been taken. All results correspond to central values of the input parameters A^L, m s (m c ), m t and ImA t . ^From 
the Table we see that the decrease in Q st corresponds to an increase in e'/e of between 21% and 63% (40 ± 11% for 
the central Bq, Bs values). The increase in the magnitude of e'/e is between 2 x 10~ 4 and 5 x 10~ 4 , to be compared 
to the current experimental world average (19.3 ± 2.4) x 10~ 4 . The magnitude of the increase will, of course, be even 
larger for models with larger values of B§. 

It is useful to comment in more detail on the application of the corrections discussed above to microscopic models 
such as the xQM and the extended NJL model [7]. Such models allow one, in principle, to compute the corrections 
corresponding to the NLO weak LEC's self-consistently within the model, as was, for example, done by the Trieste 
group [9] for Iim4 . As pointed out in Ref. [35] (in the context of the xQM), however, modifying the model predictions 
for e'/e obtained using the conventional value of Q st is more complicated than simply re-scaling the gluonic penguin 
contribution to take into account the new value of the factor 1 — tt st . The reason is that the conventional value of fl st 
(assumed to be the same for the CP-even and CP-odd cases) enters also the determination of the CP-even amplitude 
A 2 in the model; a change in il st in the CP-even sector would thus necessitate a re-fitting of the parameters of the 
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model. In fact, in Ref. [35] it was noted that the re-fitting of parameters necessitated by a shift in fl st would almost 
entirely compensate for the effect of the shifted value of Q s t in the model determination of e'/e. One should, however, 
bear in mind the caveat that this observation is based on the implicit assumption that U st is the same in the CP-even 
and CP-odd sectors. Although this is true for the non-WLEC contributions, there is no reason to expect it to be 
true for the NLO weak LEC contributions. In fact, since these contributions correspond to the hadronization of very 
different effective operators, it would be rather surprising to find them taking on the same values. Fortunately one 
docs not need to speculate idly on this question: in models such as the xQM it is possible to simply compute the NLO 
terms corresponding to the weak NLO LEC contributions. Having fitted the model parameters in the CP-even sector, 
one would then obtain, self-consistently, a determination of [Q st ] WLEC for both the CP-even and CP-odd cases. In 
order to make sure the determination of the strong IB correction to e'/e is under control, it is important to separately 
determine the NLO weak LEC contributions to the leakage amplitudes in the CP-even and CP-odd sectors. 

To summarize, we have presented a complete NLO calculation of the isospin-brcaking correction to the gluonic 
penguin operator contribution to e'/e, fl st . It is found that model-independent NLO one-loop and strong LEC 
contributions are of the opposite sign to the LO contribution, and numerically large, for typical hadronic scales. 
Combined with model estimates for the NLO weak LEC contributions, we find a significant reduction of Q s t as 
compared to the 'standard' value of 0.25 ± 0.08. Our final result is 

n st = 0.08 ± 0.05 (24) 

(where the uncertainties associated with model-dependence and BqStti have been added in quadrature). We recom- 
mend that this central value, together with, to be conservative, even larger errors, be employed in future estimates of 
e'/e in the Standard Model. 
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APPENDIX: NLO WEAK LEC CONTRIBUTIONS TO K° -> tttt 



The octet NLO weak LEC contributions to the 1 = and 1 = 2 CP-odd K° — ► tttt amplitudes are given by 

[MwLEC (^Wj ~ m l){ m K^l - "4*2") 

rex i c 2 ( 2B o(m d - m u )\ ( 2 -_ 2 - \ 
vMwlec = f2 y ^5 ) \ m K'h ~m v J 4 j 

where the are isospin-conserving NLO weak LEC combinations (the CP-even analogue of is discussed 
in Refs. [27,24]), and the Jr are isospin-breaking LEC combinations whose CP-even analogues are discussed in 
Rcfs. [19,20]. In the notation of Ref. [24] these are given by 

K- = [N~' r - 2N~' r ± 2N-' r + N~' r ] 

= [-2N-> r W 7 ~' r N-> r ± 2N w > r + AN- { r ± 2N£ r ] 
J 3 " = [Ns' r + W 6 ' r 2N-' r N^' r 4N- r - 8N-/ - 12N~/] 
J 4 " = [2N^ r ± 6N 6 ' r ± N-' r 2N w > r KWf/ 12N^ r ] . 



TABLE I. Model estimates of the NLO weak LEC combinations appearing in Eq. Al for the WDM and xQM as described 
in the text (with rv = m%/rn 2 K ). 





X QM (xlO" 3 ) 


WDM (xl0~ 3 ) 


K~ - r^K- 


-4.024 


-0.673 


J 3 ^7T J 4 


-3.571 


0.673 


TABLE II. The NLO non-WLEC contributions to Rq and R2 at the renormalisation scale /i. 


A 4 


D (non-WLEC) 


D (non-W LEC) 
It<2 




-0.01690 


-0.2359 


m p 


0.4203 


-0.3147 


TABLE III. 


The NLO weak counterterm (WLEC) contributions to the correction factors as estimated in 


the WDM and 


X QM. 








XQM 


WDM 


oWLEC 


-0.231 


-0.0331 


j^WLEC 


0.205 


-0.0331 



TABLE IV. The dependence of e'/e on Q 3 t in the Standard Model assuming, for illustrative purposes, the central values for 
Ajjj, m s (m c ), mt and ImAt as given in Refs. [12]. The units of e'/e are 10~ 4 . The values of e'/e corresponding to Q s t = 0.25 
are taken from Table 3 of the second of Refs. [12] and the range of values for Bq, Bg is the same as covered by that Table. 



B 6 


B s 


e'/e (n st = 0.25) 


e'/e (Q st = 0.08 ±0.05) 


1.0 


0.6 


8.4 


11.2 ±0.8 


1.0 


0.8 


7.0 


9.8 ±0.8 


1.0 


1.0 


5.5 


8.2 ±0.8 


1.3 


0.6 


12.8 


16.6 ± 1.2 


1.3 


0.8 


11.3 


15.1 ± 1.1 


1.3 


1.0 


9.9 


13.7 ± 1.1 
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(a) (b) (c) 




(d) (e) (f) 




FIG. 1. Feynman diagrams for K — ► mv up to 0(p 4 ) in the chiral expansion. Closed circles represent 0(p 2 ) strong vertices, 
open circles 0(p 4 ) strong vertices, closed boxes 0(p 2 ) weak vertices, and open boxes 0(p 4 ) weak vertices. No one-line weak 
tadpoles occur because, in the weak effective Lagrangian employed, they have already been rotated away. Figures (b) and (c) 
should be understood to represent collectively the strong dressing on all the external lines. 
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